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Abstract 

The Davey-Stewartson 1(DS1) system[9] is an integrable model in two dimensions. A quan- 
tum DS1 system with 2 colour-components in two dimensions has been formulated. This two- 
dimensional problem has been reduced to two one-dimensional many-body problems with 2 
colour-components. The solutions of the two-dimensional problem under consideration has been 
constructed from the resulting problems in one dimensions. For latters with the 5-function in- 
teractions and being solved by the Bethe ansatz, we introduce symmetrical and antisymmetrical 
Young operators of the permutation group and obtain the exact solutions for the quantum DS1 
system. The application of the solusions is discussed. 
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1, Introduction : 

The Davey-Stewartson 1 (DS1) system is an integrable model in space of two spatial and 
one temporal dimensions ((2+l)D). The quantized DS1 system with scalar fields (1 component 
or shortly 1C) can be formulated in terms of the Hamiltonian of quantum many-body problem 
in two dimensions, and some of them can be solved exactly Q Particularly, it has been shown 
in ref.Q that these 2D quantum iV-body system with lC-fields can be reduced to the solvable 
one-dimensional (ID) quantum iV-body systems with lC-fields and with two-body potentials Q. 
Thus through solving ID quantum iV-body problems with lC-fields we can get the solutions for 2D 
quantum iV-body problems with lC-fields. Here the key step is to separate the spatial variables of 
2D quantum iV-body problems with lC-fields by constructing an ansatzQ^] 

■ ■ ■ ,Zn,T]i,- ■ ■ > W) = Hi 1 - T e (&j) e (%'))^(6,-- ■ ,£,n)Y(t]i, ■ ■ ■ ,rj N ) 

where £jj = £j— £j and rjij = rji—rjj. This ansatz will be called the N-body variable-separation ansatz. 
It is well known that the variable-separation methods are widely used in solving high-dimensional 
one-particle problems. For instance, for geting the wave functions of electron in hydrgen atom, 
the ansatz \P(r, 9, 4>) = R(r)P(Q)§>((f)) is used to reduce the 3D problem to ID's (this ansatz is 
what we call the 1-body variable-separation ansatz). The N-body variable-separation ansatz can 
be thought of as the extension of 1-body variable-separation ansatz. Since the N-body problems 
are much more complicated than the 1-body problems, it will be highlly nontrivial to construct a 
N-body variable-separation ansatz. Ref.Q provided the first example for it and showed that the 
idea of variable- separating works indeed for the N-body problems induced from the DS1 system . 

In this paper, we intend to generalize the above idea to multicomponents DS1 system, namely 
to construct a new N-body variable-separation ansatz for the multicomponents case and to solve a 
specific model of 2D quantum DS1 system with multicomponents. 

ID N-body model with 2-components has been investigated for long timeQ|H]. The most 
famous one is the model with delta- function interaction between 2C-fermionsQ. It was solved by 
the Bethe ansatz || and leads to the Yang-Baxter equation and its thermodynamics studies |7|]|£| 
because of the completeness of the Bethe ansatz solutions. In this paper, for definiteness, we 
shall study specific 2D quantum iV-body system with 2C-fields associated with the DS1 system. 
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This quantum iV-body problem under consideration can be reduced to two ID quantum iV-body 
problems with 2C-fields of ref.Q and then be exactly solved by using an appropriate N-body 
variable-separation ansatz and the Bethe ansatz. 



2, Quantum DS1 System with Two Components in 2-Dimension : 

Following usual DS1 equation the equation for the DS1 system with two components 

reads 

^q = -^(dl + dy)q + iA 1 q + iA 2 ci, (1) 
where q has two colour components, 

-(£)■ (2 > 

and 

ip x - d y )A l = -ic{d x + 5 y )(q t q) 

(d x + d y )A 2 = ic(d x - 5 y )(q t q) 

where notation f means the hermitian transposition, and c is the coupling constant. Introducing 
the coordinates £ = x + y, rj = x — y, we have 

At = -icd^itfcd-imit) (3) 

A 2 = icfyd^tfqj + imiv) (4) 

where 

^ X (q + q) = 5 ( / H- dr ? ')q t (£,V,t)q(e, v',t), (5) 

^ J— oo Jrj 

and iti and u 2 are constants of integration. According to ref.||], we choose them as 

MO = -/W^-Oq^W-^eW,*) (6) 

MV) = IJ d^drj'U 2 ( V - V ')qH0,v',tH0,r]',t). (7) 
Thus eq.(l) can be written as 

iq = -(dl + d, 2 )q + c^a-^qtq) + fyd^tfq)]* 

+\ J dfdr/mz - O + +u 2 ( v - V)](q +, q')q, (8) 



where q' = q(£' ,7/ ,t). We quantize the system with the canonical commutation relations 

[&(e. V, t), qt(?, r/, t)} ± = 25 a M ~ ~ rf), (9) 

[qa(t,V,t)M?M,t)]± = 0- (10) 

where a, b = 1 or 2, [, ] + and [,]_ are anticommutator and commutator respectively. Then eq.(8) 
can be written in the form 

q = i[H,q] (11) 
where H is the Hamiltonian of the system 

H=\Jdidrj ( -qt(^ + ^)q + |qt[(^a-i + ^i)(qtq)] q 

+ \ J dt'drj'cfiUxiti - O + U 2 (r, - r/)]( q 't q ')q ) • (12) 
The iV-particle eigenvalue problem is 

H\V) = E\ V) (13) 

where 

|*)= I d£ l dr h ...dZ N dr }N Y, *ai...a w (6m ■ • • ZNVN)^ fam) ■ ■ ■ <&„(ZnVn) I 0). (14) 

ai...ajv 

The TV-particle wave function \& ai ... ajv is defined by eq.(14), which satisfies the TV-body Schrodinger 
equation 



Y(dl + dl)* ai ... aN + cYH^S'ivij) + efay)«'(&)]*«i. 



.ajv 



«<7 



+ E[ [/ ife) + / '^'/, ; )>I'»,..^ = S* ai ... ON (15) 

where &j = & - £j, <5'(fy ) = and e(fy) = 1 for £y > 0, for £y = 0, -1 for ^ < 0. Since 

there are products of distributions in eq.(15), an appropriate regularezation for avoiding uncertainty 
is necessary. This issue has been discussed in ref. [pi . 
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3, Variable Separation of Quantum DS1 with Two Components and Bethe Ansatz : 

Our purpose is to solve the iV-body Schrodinger equation (15). The results in ref.|| remind 
us that we can make the following ansatz 

= E IE 1 - i<tv)<vn)) 

a\ ... a N 
b[ ... b' N 

•A/" ai ... ajv ,6' 1 ...6' Jvr ^a' 1 ...a' w (^l ■ • • ^v)^...^ (Vl ■ ■ ■ W) ( 16 ) 

where M. and Af are matrices being independent of £ and r], and both X ai ... ajv (£1 . . . £w) and 
Ybi...b N (vi ■ ■ ■ Vn) are one-dimensional wave functions of N-bodies. Substituting eq.(16) into eq.(15), 
we abtian 

~^2^Ii X ai-a N +^Ui(iij)X ai ,„ aN = EiX ai ... aN (17) 

i i<j 

-E^ i y f 1 ...6 JV +E^(%)n i ... 6jv = E 2 Y bl ... bN (18) 

i i<j 

where U%(£ij) and U<i(j}ij) are two-body potentials, eqs. (17) (18) are one-dimensional iV-body 
Schrodinger equations and E\ + E 2 = E. Above derivation indicates that the two-dimensional 
TV-body Schrodinger equation (15) has been reduced into two one-dimentional TV-body Schrodinger 
equations. Namely, the variables in the two-dimentional iV-body wave function VP ai .„ OJV have been 
separated. 

At this stage A4 and J\f are unknown temporarily. It is expected that for any given pair 
of exactly solvable ID N-body problems and the correspondent solutions, we could construct the 
solutions ^Ai...A N for 2D N-body problems eq.(15) through constructing an appropriate A4 x J\f- 
matrix. It has been known that the ID N-body problem in the form of (17) or (18) can be solved 
exactly for a class of potentialsQ[^][ll]. To illustrate the construction of A4 x A/"-matrix, we take 
both potentials in (17) and (18) the delta functions Uifej) = 2gS(£ij) and UiijHj) = 2g5(r]ij) 
(g > 0, the coupling constant). Then eqs. (17) and (18) become 

-J29lX ai ... aN + 2gJ2S^ij)X ai ...a N = E t X ai ... aN (19) 

i i<j 

-J2 d l Y bi---bN+ 2 9j2 d M Y ^--bN = E 2 Y bl ... bN . (20) 

i i<j 
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As X and Y are wave functions of Fermions with two components, denoted by X F and Y , the 
problem has been solved by Yang long ago||] (more explicitly, see ref.[12] and ref.[[b]])- According 



to the Bethe ansatz, the continual solution of eq.(9) in the region of < < £q 2 < . . . < £q n < L 
reads 



X F = otp exp{i[kp 1 £ t Q 1 + . . . + kp N ^Q N ]} 
p 

12. ..N 1 21. ..If 



(Q) 
p ' 

(Q) J(kiiQ 1 +k 2 ^Q 2 +-+k N ^Q N ) , (Q) J(k2^Q 1 +k 1 ^Q 2 +...+k N ^Q N ) 
12... N ~ 21. ..N 

+(N\ - 2) others terms (21) 

where X F G {X^ ,„ „}, P = [Pi, P2, ■■, Pn] and Q = [Qi, Q2, ■ ■ ■ , Qn] are two permutations of 
the integers 1,2, . . . , N, and 

«!$... = Yjr^ Q l, (22) 
yim _ —i(kj — kj)P lm + g ^23^) 



The eigenvalue is given by 

Ex = k\ + kl + ... + k 2 N , (24) 
where {ki} are determined by the Bethe ansatz equations, 



Jkj L _ TT Kkj ~ &p) - 9/2 



-p- i(fcj - Aq) - g/2 = _ tt ijAa ~ A/a) + g , 2fi x 
l = l%-A a ) + 5 /2 l = li(A a -A^)- 5 1 j 

with a = 1, . . . , M,j = 1, . . . , N. Through exactly same procedures we can get the solution Y F and 
E 2 to eq.(20). 

As X and Y are Boson's wave- functions, denoted by X B and Y B , it is easy to be shown that 

X B = £/# ) exp{^P 1 £Q 1 +--- + «Q J v]} (27) 
p 

j(Q) _ ylma(Q) 



= ZjT^l., (28) 

7 j m ^(fcj — kj)P lm + 

Jl z'(fej - h) - g 



r?lm _ ~ kj)P lm + g 
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and the Bethe ansatz equations circ cis following |1 

ikjL _ (_]\N+i tt % ~ h + iff yr Ap — kj + ig/2 
1 ' tJikj-h-ig UAp-kj-ig/2 

TT A/3 ~ Aq + iff _ / -.xM+l A A /3 - kj + ig/2 

^Ap-Aa-ig { ' j-J^ Ap - kj - ig/2' 



(30) 
(31) 



Y B is same as X B . It is well known that X F and Y F (X B and Y B ) are antisymmetrical (symme- 
trecal) as the coordinates and the colour-indeces of the particles interchanges each other simulta- 
neously, instead of the coordinates interchanges each other merely. 

4, Young Operator of Permutation Group : 

For permutation group Sn ■ i = 1, • ■ ■ , N\}, the totally symmetrical Young operator is 



On = 5Z ei ' 

i=l 

and the totally antisymmetrical Young operator is 



^ = E(- 1 ) Ple - 



i=l 



The Young diagram for On is 



have 



1 


2 


3 




N 



N , and for An, it is 



(32) 



(33) 



. To S3, for example, we 



N 



Q _ _|_ pl2 _|_ pl3 _|_ p23 _|_ pl2 p23 _|_ p23 p!2 



_ pl2 pl3 p23 _|_ pl2 p23 _|_ p23 p!2 



(34) 
(35) 



Lemma 1: (Ojv-Xf)(£i> £2> " " , £iv) is antisymmetrical with respect to the coordinate's interchanges 
of (& < — ► £,•). 

Proof: From the definition of Oat (eq.(32)), we have 



Q N P ab = P ah O N = On- 



(36) 
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To JV = 3 case, for example, the direct calculations show 3 P 12 = P 12 = 3 , 3 P 23 = P 23 = 
3 and so on. Using eqs.(36) and (23), we have 



N Yp = {-1)O n . (37) 

From eqs.(21) and (23), X F can be written as 

X F = | e *( fe ifoi+ fc 2?o 2 +...+fejv€o^) ^ y i i2 e i(fc2SQ 1 +feiSo 2 +...+*;jvfQj V .) 

+Y^Y^e i ^i +k ^2+ k ^Q- i +-+ k N(Q N ) + (_/V! - 3) other terms N (38) 

Using eqs.(37) and (38), we obtain 

+ e i(fc2€Q 1 +fe3?Q 2 +fcl«Q 3 + ...+feiV«Q JV ) + ^1 _ 3) Qther term g }0 N a^ N 

= £(-l) p exp{i[fc ft £ 0l + • • • + kp N Z QN }}(0 N aW N ). (39) 
p 

Therefore we conclude that (OnX f )(^i, ■ ■ ■ , £at) is antisymmetrical with respect to (£j < — ► 
Lemma 2: (AnX b )(£i, • • • , £n) is antisymmetrical with respect to the coordinate's interchanges 

of (& < — > 

Proof: Noting (see eqs.(33) (29) (27)) 

A N P ab = P ab A = -A N , (40) 
A N Zfr = (-1)A N , (41) 

we then have 

(A N X B )(^,---,^ N ) = ^(- 1 ) Pe Mi[kp^Q 1 + ... + kp N ^ Q J}(A N ^ N ). (42) 

p 

Then the Lemma is proved. 

5, The Solutions of the Problem : 

The ansatz of eq.(16) can be compactly written as 

* = IK 1 " le&MVijMMXWY) (43) 

i<j 
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where {MX) and {MY) are required to be antisymmetrical under the interchanges of the coodinate 

vairables. According to Lemmas 1 and 2, we see that 

. . . r I On for ID Fermion , . lS 

M > M ={a N for ID Boson. ^ 

As the DS1 fields q a {(,v) m ec L-(l) are (2+l)D Bose fields, the commutators ([,]-, see (9) and 

(10)) are used to quantized the system and the 2D N-body wave functions denoted in ty B must be 

symmetrical under the colour-interchang {m < — > aj) and the coordinate-interchange {{^iiji) < — > 

i^jVi))- Namly, the 2D Bose wave functions ^ B must satisfy that 

>"•"•*" i.„ (45) 
As q a are (2+l)D Fermi fields, the anticommutators should be used, and '3/^ must be antisymmet- 
rical under (aj < — > aj) and ((Ci^i) < — ► (CjVi))- Namly, 

/"' "'* / ' i (46) 
Thus for the 2D Boson case, two solutions of *f? B can be constructed as following 

*? = IK 1 " z<tv)<*i))P"X F fa ■ ■ ■ M\ [o*Y F (m ■ ■ ■ w)\, (47) 

i<j 

= Hi 1 - l<&i)<riij))[MX B (t;i ■ ~tN)][AvY B ( m ■ --vn)). (48) 

i<j 

Using eqs.(36),(39),(40) and (42), we can check eq.(45) diractly. In addition, from the Bethe ansatz 
equations (25) (26) (30) (31) and E = E\ + E%, we can see that the eigenvalues of Vl/f and ^ B are 
different each other generally, i.e., the states corresponding to ^ B and ^ B are non-degenerate. 
For the 2D Fermion case, the desired results are 

*f = IK 1 " {<&j)<mj))pMX F {b ■ ■ • £*)] [A^Y B { m ■ ■ ■ tin)], (49) 

i<j 

^ = IK 1 ~ l<£vHvij))lAxX B (£i ■ ■ ■ tN)][o N Y F { m ■ ■ ■ w)]- (so) 

i<j 

Eq.(46) can also be checked diractly. The eigenvalues corresponding to are also determined by 
the Bethe equations and E = E\ + E2 ■ 

It is similar to ref.Q that we can prove an d shown in above are of the exact 
solutions of the eq.(15). Thus we conclude that the 2D quantum many-body problem induced from 
the quantum DS1 system with 2-component has been solved exactly. 
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6, The Ground-State Energies of the System : 

In this section, we discuss the ground-state energies of the DS1 system solved in the previous 
section by using the Bethe ansatz equations (25), (26) and (30), (31). Let the 2D N-body problem 
reduced from 2D DS1 system with 2 colour (or spin) components has M colours down and N — M 
colours up. Therefore both X F ' B (^, • • • &v) and Y F ' B (rn, rj 2 , ■ ■ ■ rj N ) in eqs(47)-(50) are one 
dimensional N— body wave functions with M colours down and N—M colours up. We are interested 
in the limit that N, M and the length L of the box go to infinity proportionately, i.e., both N/L = D 
and M/L = D m are finite. 

For one dimensional iV-fermion problem, by the nested Bethe ansatz (or Bethe- Yang ansatz) 

equations (25) and (26), the corresponding integration equations for the ground state readQ 

2gai(A')dA' fQi 4g Pl (k)dk 
1 " Lb i9 2 + (A-A>) 2 + L Ql g*+4(k-Ar ( > 

27TP1 = l + l- Bl g* + 4(k-Ar (52) 
where pi(k) is particle (i.e., ID fermion) density distribution function of k, and 01(A) is colour-down 
particle density distribution function of A. Namely, we have 

D= f Ql pi(k)dk, D m = f Bl <ri(A)dA, E 1 /N = D' 1 f Ql k 2 Pl (k)dk. (53) 

J-Qi J-Bx J-Qx 

For ID N-boson case, starting from the nested Bethe ansatz equations (30) and (31), similar 
integration equations for ground state of bosons can be derived (see Appendix). The results are as 
follows 

2go 2 {A')dA> [Q* 4gp 2 (k)dk 

2 ™ 2 = U 9 ^{A-A>f-Uf + ^-Ar (54) 



27TP: 



2 



_ [** Aga 2 (A)dA [Q* 2gp 2 {k')dk> 
Lb* g 2 + 4(A - kY + y_ Qa a 2 + (k — k') 2 ' lD ° j 



-B 2 g*+A{A-ky J-Q a g* + (k-kf) 2 ' 
where p 2 {k) and a 2 (A) are bosonic particle density distribution function of k and its colour-down 
particle density distribution function of A respectively, i.e., 

?2 cB 2 „ pQ 2 

p 2 (k)dk, D m = 
-Q 2 J-B 2 

The everage energies of the 2D DS1 ground states described by tyf, ^> 2 , and ^f F (see 

eqs(47)-(50)) are denoted by E(^f), E($? B ), E(9 F ) and E(^ F ) respectively. Then, the everage 
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D = f Q2 p 2 (k)dk, D m = a 2 (A)dA, E 2 /N = D' 1 [ * k 2 p 2 (k)dk. (56) 

J-Q 2 J-B 2 J-Q 2 



energies per particle for the ground-states are follows 



E(9$)/N 



E(Vf)/N 




k 2 pi(k)dk, 



p 2 (k)dk, 



(57) 



(58) 



E(9f)/N 



^(Ei + E 2 ) 
rQi 

D~ l { k 2 Pl {k)dk + 




p 2 (k)dk) 



Ei^DlN. 



(59) 



(60) 



From these equations, the follows can been seen: 1, The everage energies per particle for the ground 
states of this 2-dimensional (2D-) DS1 problem are reduced into the everage energies per particle 
of 1-dimensional (lD-)many body problems. As D and D m are given, by solving the integration 
equations (51) — (56), we obtain the pi(k) and P2(k), and then get the desired results of E(^f)/N, 
E(V%)/N, E(V{)/N and E(V$)/N. 2, For the two bosonic solutions of the 2D DS1 system 
with 2 colours (eqs (47) (48)), the everage ground state energies per particle are twice as large as 
one of lD-fermions or ID-bosons; 3, For the fermion solutions of this 2D-DS1 system, E(^>f)/N 
and E(^)/N are sum of lD-fermion everage energy per particle and ID-boson's. 4, In general, 
E($!±) / E^f) 7^ E(^i or 2 ). Namely, for same DS1 system, if the statistics of the wave functions 
(or particles) is different, the corresponding ground-state energies are different. This is remarkable 
and reflects the statistical effects in the 2D DS1 system. 

7, Discussions and Summary : 

Finally, we would like to speculate some further applications of the results presented in this 
paper to the mathematical physics. Our results may be useful in the following two respects. Firstly, 
the Bethe ansatz equations (25), (26) for fermion wave functions and (30), (31) for boson's can be 
solved respectively, even though the equations are systems of transcendental equations for which the 
roots are not easy to locate. The so-called string hypothesis is used for the analysis and classification 
of the roots for the Bethe ansatz equations QHl- Thus, we could study thier ground state, the 
excitation and the thermodynamics based on itR[@]- Then, the thermodynamical properties of 
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the ID Bose or Fermi gas with S— function interaction and with two components can be explored. 
The eqs.(47) — (50) indicate that under the thermodynamical limit the 2D DS1 gases (with two 
colour components) are classified into 2D Bose gases and 2D Fermi gases. By eqs (47) (48), the 2D 
Bose gases are composeted of two ID Fermi gases or ID Bose gases, and by eqs (49) (50), the 2D 
Fermi gases are composeted of ID Fermi gas and ID Bose gas. Hence, the thermodynamics of 2D 
DS1 gases with two colour components can be derived exactly. It would be interesting in physics, 
because this is an interesting and nontrivial example to illustrate coupling (or fusing) of two ID 
2-component gases with 5— function interactin and with different or same statistics. Secondly, the 
colourless DS1 equation originated in studies of nonlinear phenomena S. Five years ago, Pang, 



Pu and Zhao [14] showed an example that the solutions of the initial-boundary- value problem for 
the related classical DS1 equation in ref.|l5| are consistent with the solutions for the quantum 
DS1 system with time-dependent applied forces. This indicates that the classical solutions of DS1 
equation are corresponding to the classcal limit of the solutions for the quantum DS1 system. This 
is actully a new method to reveal the solutions of the colourless DS1 equation. To the quantum DS1 
system with colour indices studied in this present paper, similar correspondences are expectable. 
Hence, the structure of the solutions of the quantum DSl-system with colour indices revealed in 
this paper would be helpful to understand the corresponding classical solutions of DS1 systems 
with colour. The specific studies on the above speculations would be meaningful, however, they 
are beyond the scope of this present paper. 

To summarize. We formulated the quantum multicomponent DS1 system in terms of the 
quantum multicomponent many-body Hamiltonain in 2D space. Then we reduced this 2D Hamil- 
tonain to two ID multicomponent many-body problems. As the potential between two particles 
with two components in one dimension is 5— function, the Bethe ansatz was used to solve these 
ID problems. By using the ansatz of ref.[jl| and introducing some useful Young operators, we pre- 
sented a new N-body variable-separation ansatz for fusing two ID-solutions to construct 2D wave 
functions of the quantum many-body problem which is induced from the quantum 2-component 
DS1 system. There are two types of wave functions: Boson's and Fermion's. Both of them satisfy 
the 2D many-body Schrodinger equation of the DS1 system exactly. The results have been used to 
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study the ground states of the system. Some further applications of the results presented in this 
paper are speculated and discussed. 
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Appendix 

Let us derive eqs (54) and (55) in the text. We start from the Bethe ansatz equations (30) and 
(31) of ID bosons with two colour componts. Taking the logarithm of (30) and (31) respectively, 
we have 

N , , M 



kjL = 2tt4 - 2 £ tan" 1 ^— ^ - 2 £ tan" 1 2(A/3 k]) 



i=i 9 — x 9 



(Al) 



2 £ tan" 1 A/3 A " = 2vr J A + 2 £ tan" 1 2(A/3 kj \ (A2) 



a=l ^ ^ 

where (for the case of A =even, M =odd) 

^ + 7fe = successive integers from 1 — ^A to + ^A, 

Ja = successive integers from — -^(M — 1) to + -^(M — 1). 

We can now approach the limit A — ► oo, M — > oo, L — > oo proportionally, obtaining 

k = 2tt/ 2 - 2 [ Q2 dk'p 2 (k') tan" 1 ^ - k>) - 2 dAa 2 (A) tan" 1 2 ( A ~ fc \ (A3) 
J-Q 2 9 J-b 2 g 

rB 2 A — A' rQ2 of A — frl 

2 / dAV 2 (A') tan" 1 = 2nh 2 + 2 / dfc^fc) tan" 1 -i ^, (A4) 

•/-B 2 5 J-Q2 9 

dh 2 df 2 

-dA =<72 > ~dk=^ {Ab) 

A /"Q2 ,,,„ „ M f B * 
p 2 {k)dk, D m = — = \ 
-Qi L> J-B 2 

Or, after differentiation, 



A f Q2 M f B2 

D = -= P2 (k)dk, D m = — = a 2 (A)dA. (A6) 
L J-n 2 L J-Bo 



_ r B * 2ga 2 (A')dA' [Q> Agp 2 {k)dk 
a2 ~ J-B 2 g 2 + (A- A') 2 J. Q2 g 2 + 4(k-Ar ( ?) 

J-b 2 g 2 + 4(A - k) 2 J-Q 2 g 2 + (k-k>) 2 ' 
which are just eqs (54) and (55). 
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